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Abstract. We have made cross-correlation analyses of (2 
- 15 keV) HEAO-1 A2 and 1 keV ROSAT PSPC All- 
Sky Survey maps over a selected area (~ 4000 deg 2 ) with 
high galactic latitude (6 ^ 40°). We have calculated the 
correlations for the bright ROSAT sources and residual 
background separately with the HEAO-1 A2 TOT (2 - 10 
keV) and HRD (5-15 keV) maps. In any case, the angular 
dependence of the cross-correlation function (CCF) was 
consistent with expected from the HEAO-1 A2 collima- 
tor response function and the smoothing function of the 
ROSAT 1 keV map. The amplitude of the bright ROSAT 
source - A2 CCFs are consistent with expectations from 
model populations of AGNs and clusters of galaxies, which 
emit in both bands. However, the residual ROSAT back- 
ground - A2 CCFs amplitude at zero degree are about a 
factor of three larger than that expected from the model 
populations. 

Our soft-hard zero-lag and angular CCF results have 
been compared with the 1 keV auto-correlation function 
(ACF) found by Soltan et al. (1995) for the same ROSAT 
data. Their significant angular CCF at a scale of 9 55 5° is 
of an unknown origin and it could represent a new class of 
cosmic X-ray sources. If this 1 keV ACF has a hot plasma 
spectrum with kT ~ 2 keV, contribution of this compo- 
nent is consistent with both our zero-lag CCF in excess of 
the population synthesis model prediction and the upper- 
limit to the angular CCF at 9 ~ 2.5°. On the other hand, 
if this component has a lower temperature or a steeper 
spectrum, a major modification to the population synthe- 
sis model and/or an introduction of new classes would be 
needed. 
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1. Introduction 

With deep ROSAT observations of blank fields, up to ~ 
60% of the soft (^ 2keV) cosmic X-ray background (XRB) 
has been resolved into individual sources (e.g. Hasinger 
et al. 1993). A number of groups have worked/are cur- 
rently working on identification programs of these ROSAT 
sources and we can expect that we will have the over- 
all picture of the origin of most of the soft (~ 1 keV) 
X-ray background including contributing source popula- 
tions and their cosmological evolution. With ASCA ob- 
servations of blank fields, the XRB at higher energies (2 - 
10 keV) have also been resolved into sources, although the 
resolved fraction of the total XRB flux is smaller than that 
of ROSAT because of the limited spatial resolution. While 
these observational efforts are rapidly in progress, the uni- 
fied theory of AGNs and supporting X-ray observations of 
AGNs for the basic framework of the unified scheme have 
provided insights to the origin of the bulk properties of 
the XRB. In particular, intrinsic photoelectric absorption 
observed in the X-ray spectra of Seyfert 2 galaxies (e.g. 
Awaki et al. 1991; Mulchaey et al. 1992), consistent with 
the view of the unified scheme that a Seyfert 2 galaxy is a 
Seyfert 1 seen at an angle where a torus surrounding the 
central engine blocks the line of sight. A wide range of ab- 
sorption column densities have been observed (Schartel et 
al. 1995) for a hard X-ray flux-limited sample of AGNs by 
Piccinotti et al. (1982). Models of the XRB with evolving 
populations of unabsorbed and internally-absorbed AGNs 
with various absorption column densities, based on the 
unified scheme of AGNs have been successfully made con- 
sistent with many observational properties of the XRB 
(Madau et al. 1994; Comastri et al. 1995). In particular, 
they explained the global spectral shape of the XRB, es- 
pecially the characteristic break at 30 keV and source 
counts in soft and hard X-rays. While these models are 
basically successful, there is still room for major modifica- 
tions, given current observations. A recent ASCA observa- 
tion of NGC 3628 by Yaqoob et al. (1995) showed that the 
X-ray (E 10 keV) spectrum of this spiral galaxy, which 
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had previously been classified as a starburst, is rather flat. 
Based on this observation, they argued that the contri- 
bution of low-luminosity extragalactic X-ray sources, like 
NGC 3628, to the E ^ 2 keV XRB could be significant. 
Such a low X-ray luminosity population has been deliber- 
ately neglected in the Madau et al. (1994) and Comastri 
et al. (1995) models. Also these models assumed a sim- 
ple version of the AGN unified scheme where the ratio 
of the numbers of absorbed and unabsorbed objects does 
not depend on their intrinsic X-ray luminosity. However, 
Barcons et al. (1995) argued that this is not likely to be 
the case but rather the fraction of the unabsorbed AGNs 
(type 1) should increase with the intrinsic luminosity. 

Besides detailed studies of the nature of individual 
X-ray sources which can contribute to the XRB, sta- 
tistical properties of the spatial fluctuation of the un- 
resolved background provide us with other perspectives 
of cosmic X-ray sources. Fluctuations of the unresolved 
XRB have been used to infer the log TV — log S rela- 
tion below the source detection flux-limit (e.g. Shafer 
1983; Hasinger et al. 1993; Barcons et al. 1994). The an- 
gular auto-correlation function is often used to describe 
spatial structure of the XRB. Many authors discussed 
the constraints on the combination of the X-ray volume 
emissivity and clustering properties of the X-ray sources 
from the ACF (or upper limits of ACF) (e.g. Danese 
et al. 1993; Carrera et al. 1991; Soltan & Hasinger 1994). 
While the XRB ACF traces the properties of all the X- 
ray sources along the line of sight and thus the interpreta- 
tion is highly model dependent, cross-correlating catalogs 
of known sources with the XRB gives more concrete in- 
formation on the XRB constituents. In particular, cross- 
correlation functions (CCF) of unresolved XRB spatial 
fluctuations with galaxy catalogs give information on the 
contribution of these galaxies and objects clustered with 
these sources to the XRB (Jahoda et al. 1991, 1992; Lahav 
et al. 1993; Miyaji et al. 1994; Roche et al. 1994; Barcons 
et al. 1995). 

Auto and cross-correlation function analyses over a 
large region of the sky occasionally discover new compo- 
nents of cosmic X-ray emission. The auto-correlation of 
the HEAO-1 A2 hard X-ray map nearly over the whole sky 
at high galactic latitudes revealed a weak large-scale com- 
ponent at ^ 40°, which perhaps is associated with struc- 
tures in the supergalactic plane (Jahoda 1993). Using the 
ROSAT All-Sky Survey (RASS) data, Soltan et al. (1995) 
found an extragalactic extended component in the ACF 
of at 1 keV at a scale of 5 degrees. They also found an 
angular correlation at 5 - 10 degree scales in the CCF 
between the 1 keV XRB and the positions of the Abcll 
clusters. They modeled the latter component as a 10 Mpc 
scale diffuse gas surrounding clusters of galaxies, which it- 
self had not been recognized in any previous observation. 
They found, however, that clusters of galaxies, consider- 
ing this putative circum-cluster emission and clustering of 
clusters, explain only about a third of their extended 1 keV 



ACF component. Thus the main portion of the extended 
ACF could represent still a new population of cosmic X- 
ray emission and investigating their nature is only possible 
by statistical analyses over a large region of the sky. 

In this paper, we present the results and possible in- 
terpretations of a cross-correlation analysis between the 
soft (-1 keV) X-ray sky observed with ROSAT PSPC 
All-Sky Survey and the hard (2 - 15 keV) X-ray sky ob- 
served with the HEAO-1 A2 experiment. The purpose of 
the correlation between soft-hard X-rays presented in this 
paper is two-folded. Zero-lag and angular cross-correlation 
functions (CCFs) between surface brightnesses contain in- 
formation on X-ray sources common to these two bands. 
Thus this would give additional constraints to the popu- 
lation synthesis models. Especially, it would give a check 
to the models based on unabsorbed and absorbed AGNs. 
A large area correlation study also may give a clue to 
the nature of the extended ACF component that Soltan 
et al. (1995) found. The scope of this paper is as follows: In 
Sect.|], we explain the data used. We explain the CCF cal- 
culation and present the results of the calculation in Sect. 
[| Sect^ presents formulations relating observed CCF 
with models including individual sources in the population 
synthesis models and extended components. In Sect.||, we 
use the formulation to compare various models with the 
observations. The purpose of the section is to check the 
effects of various components on the observed CCF and 
we do not intend to construct a complete model consis- 
tent with all the existing constraints. Such an elaborate 
modeling would be a topic of a future paper. Finally we 
conclude our discussion in Sect. |^. 

2. Data 

2.1. The ROSAT All-Sky Map 

We have used a clean surface brightness map constructed 
from the ROSAT (Triimper 1983) PSPC (Pfeffermann 
et al. 1986) All-Sky Survey (Snowden & Schmitt 1990; 
Voges 1992) in the north galactic pole region. In mak- 
ing the map, the period of Short-Term Enhancements 
has been excluded and non-cosmic backgrounds (parti- 
cle background, solar scattered X-rays and the Long- 
Term Enhancements) have been subtracted as explained 
in Snowden et al. (1995). Soltan et al. (1995) selected a rel- 
atively large portion of the high galactic latitude sky which 
is free from contamination by galactic structures to inves- 
tigate the extragalactic structures in the soft X-ray sky: 

70° < I < 250°, b > 40°. (1) 

We also use this area for our correlation with the hard X- 
ray sky. We also included a 3.5 degree wide strip surround- 
ing this region in our analysis, taking the larger HEAO-1 
A2 collimator response function into account. In this work, 
we have only used the map in the R6 (PI channel range 
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91 - 131, corresponding to 0.9 - 1.3 keV) band (Snowden 
et al. 1993), which is least contaminated by the galactic 
emission, scattered solar X-rays, and particle background. 
The energy response curve of the R6 band is shown as a 
solid line in Fig. |l|. 

The conversion factor between the observed R6 count 
rate and the 0.5 - 2 keV flux before the absorption by 
the Galaxy is 3.24xl0 -11 erg s _1 cm _2 /(R6 cts s _1 ) for a 
power-law photon index of (a p h — 2.0) for the absorption 
column density of TVh = 1.8 X 10 20 cm~ 2 , which is the 
average galactic value for the region. The range of the 
galactic column in the region is 0.6 — 4. x 10 20 cm~ 2 and 
the conversion factors in that region differ from the mean 
value by ^ 5%. Therefore we use the average Nh value for 
the whole region for the analysis. 

t — i i i | 1 1 1 — i — i — i i i | 1 1 1 — i — r 




1 10 
Energy (keV) 



Fig. 1. Energy response curves of the ROSATR6 and HEAO-1 
A2 TOT/HRD bands, used in our CCF analysis, are shown as 
labeled. Each curve is normalized at its peak 



2.2. The HEAO-1 A2 Hard X-ray Map 

We have used the all-sky hard X-ray map constructed 
from the combination of the HED3 and MED detectors 
of the HEAO-1 A2 experiment (Rothschild et al. 1979). 
The combination of the detectors is sensitive in 2-60 keV. 
Conventionally all-sky maps from this combination are an- 
alyzed in a number of standard colors defined as weighted 
sums of counts from certain combinations of detector lay- 
ers and pulse height channels (discovery scalers; see Allen 
et al. 1994). In this work, we have used the TOTAL (in- 
dicated by TOT) and HARD (indicated by HRD) band 
maps. The energy response curves of the TOT and HRD 
bands are shown as dashed and dotted lines respectively 
in Fig. |IJ For the cosmic X-ray background and canonical 
AGN spectrum, most of the photons come from the 2-10 
keV range for the TOT map and 5-15 keV range for the 



HRD map. For a power-law spectrum with photon index 
<3ph = 1-7, the conversion factor between the TOT count 
rates and the 2-10 keV flux is 2.2 x 10~ n erg s^ 1 cm~ 2 
(TOT cts s -1 ) -1 and the conversion factor between the 
HRD count rates and the 5-15 keV flux is 1.8 x 10~ n 
erg s _1 cm~ 2 (HRD cts s -1 ) -1 . These conversion factors 
vary only by 4% for a 40 keV bremsstrahlung spectrum, 
an a = 0.7 power-law, the same power-law with neutral 
gas absorption up to Nn ^ 10 24 cm~ 2 . 

The maps used here have the beam profile (or the 
collimator response function) well represented by (Shafer 
1983): 

5 A2 =nmx(l-M,0)max(l-M,0) (2) 

where y represents the coordinate along the scan direction 
of the survey, which is along the great circle of a constant 
ecliptic longitude, and x is the coordinate along the direc- 
tion perpendicular to the scan path. 

2.3. Excluded Regions 

Regions around conspicuous sources in the R6 band have 
been excluded from the analysis in order to avoid the sit- 
uation that a few sources dominate the correlation signal. 
Those are the Coma cluster, Abell 1367, Mkn 421 and 
NGC 4253 (Mkn 766). These are all the sources in the re- 
gion with R6 count rates greater than 0.9 counts per sec- 
ond. We have excluded 3? 4 radius regions (corresponding 
to the maximum extent of the A2 beam) around Mkn 421 
and NGC 4253. The exclusion radius around the two con- 
spicuous galaxy clusters (i.e. the Coma cluster and Abell 
1367) was 7°. In addition, we have also excluded the 3.4 
degree radius regions around the HEAO-1 A2 sources in 
Piccinotti et al. (1982). The flux cutoff of this complete 
flux-limited catalog is 1.25 counts s _1 in the HEAO-1 A2 
R15 band, which is similar to the TOTAL band used here 
but a little harder (Allen et al. 1994). 

3. The Soft-Hard Cross- Correlation Function 

3.1. The CCF calculation 

We have calculated the angular cross-correlation functions 
between the ROSAT All-Sky Survey R6 and HEAO-1 A2 
maps in the region defined in Eq. ([!]). Since the ROSAT 
All-Sky map has a much higher spatial resolution and 
denser samplings, we have smoothed the ROSAT data be- 
fore correlating. As the smoothing function, we took the 
form represented by Eq. (||) but smaller FWHMs in or- 
der to improve the signal-to-noise ratio of the correlation. 
Among a number of sizes we tried for the smoothing func- 
tion, the best results could be obtained when we chose 1?5 
and 0?75 (FWHMs) in the x and y directions respectively. 

As a statistical measure to indicate the degree of cor- 
relation between two maps (here we denote the two maps 
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by subscripts s and h, representing soft and hard respec- 
tively), we calculate the angular cross-correlation function 
(CCF): 



W sh {0) 



(Is) (h) 



(3) 



where I s and I h are the ROSAT and HEAO-1 A2 inten- 
sities measured at the gridding points which are within 
the region defined above. We have taken pairs along the 
scan path (ecliptic longitude) to calculate the angular cor- 
relation in order to make use of the narrower collimator 
response. 

The uncertainties of the correlation functions are esti- 
mated by correlating the ROSAT map with different parts 
of the sky in the HEAO-1 A2 map. We made this by 
rotating the ROSAT map around the galactic pole and 
calculating the correlation with the A2 map in the same 
way. We have also transposed the ROSAT map about the 
galactic equator followed by rotations about the galactic 
pole to calculate the correlation with the A2 map. Thus 
we produced a number of different artificial null samples 
with equivalent statistical properties. Regions around Pic- 
cinotti et al. sources are excluded as we did in the real cor- 
relation calculation. We also excluded 20° and 25° regions 
around LMC and SMC respectively. We have rotated the 
map in steps of 30° on both hemispheres, thus we get 23 
null-hypothesis CCFs between unrelated parts of the sky. 
These give a fair estimate of the uncertainties of the cor- 
relation function. 

3.2. CCF of ROSAT bright sources with the HEAO-1 A2 
Maps 

In order to evaluate the contribution to the CCF from 
bright ROSAT sources, we have made a list of bright 
sources by performing a source detection procedure to the 
R6 map in the region defined in Eq.(0). The source de- 
tection has been made to the map gridded in 12' x 12' 
pixels, larger than the ROSAT PSPC point spread func- 
tion, and thus is much less sensitive than the source de- 
tection procedure directly applied to the full-resolution 
ROSAT All-Sky Survey data. A ROSAT All-Sky Survey 
bright source list using such a procedure (Voges 1995) is 
in progress at the time of writing this paper. An analysis 
using the source list will be a topic of a future paper. In 
this work, we use the coarsely rebinned map to obtain a 
complete R6 flux-limited sample of bright sources. As de- 
scribed by Hasinger et al. (1993), the source detection have 
been made in two steps: the first step is to slide a detection 
cell over the map and find excess over the region imme- 
diately surrounding the cell (LDETECT). A background 
map was created by removing the sources detected in the 
LDETECT procedure, filling the holes by the mean counts 
surrounding them, and making a bi-cubic spline fit to the 
filled map. This procedure gives a more robust measure 
of the backgound for the source detection than taking the 



background from immediate surroundings. Then we slid 
the detection cell again to the original map to find signif- 
icant excess counts over the background map created in 
the above process (MDETECT). 

The source detection have been made using procedures 
within the EXSAS software package (Zimmermann et al. 
1994). The size of the detection cell has been taken to be 
3x3 cells (36' x 36') and the source detection threshold 
is set at a significance of — log Pi = 10, where Pi is the 
probability that the random Poisson fluctuation make the 
excess counts in the cell. This roughly corresponds to a 4<r 
excess of a Gaussian fluctuation. Since the i?05^4TPSPC 
PSF and extents of clusters of galaxies (except for closest 
ones, which have been excluded from the analysis (Sect. 
2.3 ) are much smaller than the size of the detection cell, 
the source list from the MDETECT procedure is expected 
to be free from significant biases. 

In the least sensitive areas, the detection limit was 0.08 
R6 cts s . If we assume a Seyfert 1-like broken power-law 
spectrum with a p h = 2.3 for E < 1.5 keV, a p h = 1.7 for 
E > 1.5 keV, and the average galactic JVh of the region 
(= 1.8 x 10 20 cm" 2 ), 0.08 R6 cts s" 1 corresponds to a 0.5 
- 2 keV flux of 2.7xl0" 12 ergem"^" 1 and 2-10 keV 
flux of 3.8xl0~ 12 erg cm" 2 s _1 . 

By this procedure, we have obtained a complete flux- 
limited (> 0.08R6 cts s" 1 ) sample of 183 sources above 
this limit over the 4575 deg 2 area of the sky (0.04 
sources deg -2 ). The log N - log S function of the de- 
tected sources has a Euclidean slope. The number den- 
sity corresponds to that of the Einstein Medium Sensitiv- 
ity Survey (EMSS) (Gioia et al. 1990) above 4.1xl0" 12 
erg cm" 2 s" 1 in 0.3 - 3.5 keV. 
the 0.3 - 3.5 keV flux versus 
sion factor is ~ 5 x 10" 11 ergs""cm 
For a power-law spectrum with a p h 



They are consistent if 
R6 count rate conver- 
i„™-2 (Rectss" 1 )" 1 . 
2.0, this factor 

is ~ 5.7 x 10" 11 ergs" 1 cm" 2 (R6 cts s" 1 )" 1 and for a 
kT = 2 keV Raymond-Smith spectrum (heavy element 
abundance=0.4 cosmic), this number is ~ 4.4 x 10" 11 



erg s cm (R6 cts s 



Thus our bright source counts 



are consistent with the EMSS counts considering a possi- 
ble range of the source spectra. 

We have calculated the CCF between the detected 
sources (R6 count rate weighte d) a nd the HEAO-1 A2 



maps in the same way as Sect. 3.1. The errors of these 



CCFs have been calculated also in the same manner. 

Using the detected sources, we have also constructed 
a source-removed ROSAT map. This has been made by 
replacing the 3x3 image pixel regions around the sources 
with R6 count rate greater than 0.08 cts s _1 by the pixel 
values in the same regions of the background map pro- 
duced for the MDETECT procedure (see above). We have 
also calculated the CCF between the source-removed map 
and the HEAO-1 A2 maps. The source-removed map still 
contains sources below the limit and we expect that such 
fainter sources still contribute to the correlation. 
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Table 1. The Zero- lag Correlation Results 



Id. 


R6 cr. range 

[cts s" 1 ] 


HEAO-1 
A2 Band 


/ 1 \ 

(L) 

[cts s _1 


/ r \ 
(Ih) 

dee" 2 l 




Wsh(0) 




Tl 


- 0.9 


TOT 


2.9 x 10" 1 


7.2 x 10" 


i 


(3.2 ±0.4) x 10" 


3 


T2 


- 0.08 


TOT 


2.9 x 10" 1 


7.2 x 10" 


i 


(1.7 ±0.3) x 10" 


3 


T3 


0.08 - 0.9 


TOT 


6.1 x 10 _3a 


7.2 x 10" 


i 


(6.6 ± 1.0) x 10" 


2 


T4 


0.2 - 0.9 


TOT 


2.6 x 10" 3a 


7.2 x 10" 


l 


(8.3 ± 1.7) x 10" 


2 


HI 


- 0.9 


HRD 


2.9 x 10" 1 


8.6 x 10" 


l 


(2.4 ±0.4) x 10" 


3 


H2 


- 0.08 


HRD 


2.9 x 10" 1 


8.6 x 10" 


l 


(1.0 ±0.3) x 10" 


3 


H3 


0.08 - 0.9 


HRD 


6.1 x 10 _3a 


8.6 x 10" 


l 


(5.4 ± 1.6) x 10" 


2 


H4 


0.2 - 0.9 


HRD 


2.6 x 10 _3a 


8.6 x 10" 


l 


(6.6 ±2.6) x 10" 


2 


a 


Intensity from 


sources in 


the quoted flux 


range 









3.3. Results of the CCF Calculations 

Significant zero-lag correlation signals have been de- 
tected in most of the cases we have calculated (9 = case 
of Eq. @). Table summarizes the zero-lag correlation 
results for the calculations we made for various cases. The 
second column of Table [l] is the R6 count rate range of 
the sources correlated with the A2 map. A zero value of 
the lower bound in this column means that the R6 sur- 
face brightness map has been used and a non-zero value 
represents that th e s ource list in the quoted flux range 
explained in Sect. B.2 has been used. 



The angular correlation functions for selected cases 
have been plotted in Fig. |^ along with the angular cor- 
relation functions between the rotated R6 maps with the 
A2 map to show the dispersion of the CCF when no real 
sky-correlation exists. The label in each figure corresponds 
to the correlation ID in Table [|. Also the angular depen- 
dence of the CCF in case of the purely Poissonian case (see 
Sect. f|) due to the A2 collimator reponse and the ROSAT 
All-Sky Survey map smoothing, normalized at 9 = 0, has 
been plotted on each plot. In any case, we did not find sig- 
nificant extended component in the angular CCF beyond 
expected from the Poisson effect. 




(a) R6 BG vs A2 TOT (T2) 
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(b) R6 SRCs vs A2 TOT (T3) 



(c) R6 BG vs A2 HRD (H2) 
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I (d) R6 SRCs vs A2 HRD (H3) 
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Fig. 2. The angular CCF for selected cases as labeled. The 
open octagons show the observed correlation and the dotted 
lines show the 23 uncorrelated cases. The solid line is the an- 
gular dependence of the Poisson effect caused by the A2 colli- 
mator response function and smoothing 



4. Basic Formulations 

The CCF depends on Poisson and clustering effects. The 
Poisson effect is caused by individual sources which emit 
both soft and hard X-rays and the angular dependence 
of the correlation coefficient simply reflects the overlap of 
the finite-sized beams of soft and hard X-ray observations 
and smoothings. The clustering effect reflects the real sky 
structure of the X-ray sources common to (or at least cor- 
related with) both soft and hard X-rays. This structure 
may be due to clustering of sources or a diffuse X-ray 
component which extends to a scale comparable or larger 
than the size of the beam. If one is concerned with real 
extended structure, the distinction between Poisson and 
clustering effects is somewhat ambiguous. In our case, the 
extent of the X-ray emission from an individual galaxy 
cluster is much smaller than the beam size of the HE A O- 
1 A2 experiment and thus we treat them in the Poisson 
term. On the other hand, we treat the effect of possible 
diffuse components with scales of several degrees in the 
clustering effect. 

The observed angular correlation function is then, ex- 
pressed by the sum of two terms (cf. Lahav et al. 1993; 
Miyaji 1994; Miyaji et al. 1994): 



W sh (9)(I s )(I h )=r )p (e)+r 1c {0), 



(4) 
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where T s ^ represent the intensities of soft and hard X-rays 
per beam respectively: 



</ s , h )B s , h (f2)df!= (7 Sjh )ft 



Bs,h, 



(5) 



where Q is the unit vector towards the solid angle clement 
dO. 

As we will see in Sect. ||, the correlation is dominated 
by local sources and cosmological evolution is not impor- 
tant in our case, because of the large beam of the HEAO-1 
A2 experiment. However, we will develop a full cosmolog- 
ical formulation below for completeness. 

4-1- Poisson Effect from Common Sources 

We can express the Poisson term of the cross-correlation 
function as: 

V P (0)= J J S s S h N sh (S s ,S h )dS s dS h x 

( [ B B (Cl)B h (Cl-0)dn), (6) 



where iV s h(S s , Sh) dS s , dSh is the bivariate N(S) function 
defined as the number of sources per solid angle whose soft 
and hard X-ray fluxes fall in the ranges [S s — 0.5 dS s , S s + 
0.5 dS s ] and [Si, -0.5 dS h , S h +0.5 dS h ] respectively. Also 9 
is the offset vector between hard and soft X-ray measure- 
ments. The flux integrations (over S s or Sh) are over source 
fluxes included in the corresponding maps. When we use a 
surface brightness map where bright sources above a cer- 
tain limit are excluded, this limit defines the upper bound 
of the integration and the lower bound is zero. When we 
use a complete flux-limited source list, this flux limit de- 
fines the lower bound. 

The double integral over the soft and hard fluxes 
in Eq. (6) can be calculated for given population syn- 
thesis models. Let us suppose we have evolving popula- 
tions of objects with the bivariate soft-hard luminosity 
function (for the i-th population) per comoving volume: 
&i(L s , Lh, z) dL s dLh, where each population is assumed 
to have a universal spectral shape for simplicity. Then, 



S s S h N(S S , S h ) dS s dS h 



16w 2 H Jo 



D L (z)2(i + z )3(i + n o2 )i ± poi 



r [K si (z)K hi (z): 



L s L h $i(L s ,L h ,z)dL h dL s } dz, 



(J) 



L m.n 



where Dl(z) is the luminosity distance as a function of 
redshift and if s (h)i is the soft (hard) band K-correction 
for the i-th population. The luminosity integrations are 
over luminosities where the observed flux falls in the range 
of source fluxes included in the corresponding correlating 



maps. Thus the bounds of the luminosity integrations de- 
pend on redshift and the spectrum of the population. Now 
let the soft and hard X-ray luminosities of each population 
be related by L s — aihh- Then the double integral over 
hard and soft X-ray luminosities can be reduced to: 



h 

min(L™ ax ,a^L 



max 



(if 



L s Lh§i(L s , Lh, z)dL\ l dL s 
di L 2 S $ s (L s ,z)dL s , 



(8) 



where <fr s (L s ,z) is the comoving soft X-ray luminosity 
function, which we can find in literature based on ROSAT 
surveys (e.g. Boyle et al. 1993, 1994). In the case of the 
pure-luminosity evolution of the power-law form, L(z) = 
L(0)(1 + z) p , the evolving luminosity function can be ex- 
pressed by the present epoch one: 



$> si (L s , z) = $ si ([l + z}~ p L s , 0)(1 + z)~p. 



.2. Extended Component 



(9) 



The underlying angular CCF of the soft and hard X-ray 
skies 



_ (Is(Cl)h(Cl + 6)) 

Wsh{V) = TTTTT^ 1, 



> 



(10) 



(Is){h) 

reflects the extended structure common to soft and hard 
X-rays and clustering of X-ray sources. 

The observed angular correlation function (W s h(0)) is 
contributed by the underlying angular cross-correlation 
function w s h(0) of the real sky. The term r] c in Eq. |] re- 
flects this component: 



Vc (9) = (I s )(I h ) x 



(11) 



where 8' = |fii — Q2 + 6\. In effect, the underlying angular 
correlation function is smoothed with the convolution of 
soft and hard beams. At large separations (9 is much larger 
than the scale size of the beam), W s h(9) ~ w s h(9) and the 
Poisson term is zero. 

The underlying w s h can be expressed by the cosmolog- 
ically evolving soft and hard volume emissivities and the 
spatial correlation function of soft and hard X-ray sources. 
The expression is quite parallel to the case of ACF found 
in literature (e.g. Danese et al. 1993; Soltan & Hasinger 
1994). Let us consider populations of soft and hard X-ray 
emitting sources with redshift dependent comoving vol- 
ume emissivities of p s (z) and f>h(z). Let the spatial cor- 
relation function between these populations be £ s h(f*> z). 
Then, in the small-angle long-distance approximation, 



Wsh (0) 



Wir 2 H 



p s (z)Ks(z)ph(z)K h (z) x 



U(V{Da9) 2 +x 2 , z) dx{l + z)- 4 (l + r> z)5dz, (12) 
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where K s (z) and K^{z) are the K-corrections for the 
soft and hard X-ray bands respectively and Da{z) (— 
Dl(z)(1 ± z)~ 2 ) is the angular distance. 



because this quantity is proportional to the contribution 
from each of the R6 sources at a certain flux-range and 
source-removed background. 



5. Comparison of Models with Observations 

5.1. AGN and Cluster Contributions 

We have calculated the expected contributions from two 
major classes of extragalactic X-ray sources, i.e. AGNs 
and clusters of galaxies, to the observed zero-lag correla- 
tion strengths using the population synthesis technique. 
For the following models, the standard cosmology with 
A = and qo = has been used. We also denote the Hub- 
ble constant by Ho — 50/i5okms _1 Mpc -1 . Note that the 
expected correlations from the models do not depend on 
the value of the Hubble constant. 

As the first model (Model PI, where P represents 
the Poisson effect) we have constructed an AGN popu- 
lation synthesis model following the recipe by Comastri 
et al. (1995) for their baseline model. Contributions from 
X-ray clusters of galaxies have been added. In adding the 
cluster contribution, we have used the analytical forms 
of the 2-10 keV X-ray luminosity function (XLF) by 
Edge et al. (1990). We have used their volume-limited 
expression for z < 0.1 and flux- limited expression for 
0.1 < z < 0.17, taking into account the deficit of most 
luminous clusters at higher redshifts they observed. In 
the 0.17 < z < 0.6, we have used the power-law forms 
of the XLF given by Henry et al. (1992) in their three red- 
shift bins. They showed the steepening of the XLF for 
higher redshift bins. Thus their power-law form would pre- 
dict too large an XLF at low X-ray luminosities (below 
their sample limit). Thus we take the smaller (in comov- 
ing coordinates) of the Edge et al. analytical form and 
the Henry et al. power-law form. We have taken the 2- 
10 keV X-ray luminosity range of 0.01 < £ X 44 < 120, 
where L X 44 represents X-ray luminosity measured in units 
of 10 44 /i^ 2 ergs _1 . The clusters have been assumed to 
have X-ray spectra represented by the Raymond & Smith 
plasma with a heavy element abundance of 0.4. The tem- 
peratures of a cluster have been chosen from kT = 1, 
1.5, 2.5, 4, 7, 10 keV, according to its luminosity between 
the minimum and the maximum values given above. The 
dividing luminosities are Z/ X 44 = 0.02, 0.1, 0.5, 3.0, and 
15.0, which roughly corresponds to the X-ray luminosity 
versus temperature correlation (David et al. 1993). There 
is a significant scatter in the luminosity - temperature 
correlation but this has not been taken into account in 
the model. However, the results are insensitive to the de- 
tailed assumption on this L x -T correlation. For example, 
assuming a single temperature of kT = 4 keV for all clus- 
ters changed the result by 10% for T2, which is somewhat 
smaller than the statistical error of the observation. The 
expected correlation from this model is compared with 
data in Table 0. We compare the values of W s h(0) (J s ) (Jh) 



Table 2. Comparison of correlation results with models 



W sh (0)(I s }{h) 
Id. [1CT 4 R6 cts s" 1 A2 cts s" 1 ] 





Observed 






Models a 






PI 


P2 


P3 


El 


E2 


T2 


3.5 ±0.6 


0.9 


1.3 


0.9 


3.1 


1.9 


T3 


2.9 ±0.4 


2.5 


3.6 


2.5 






H2 


2.5 ±0.7 


0.7 


1.1 


0.8 


0.8 


0.2 


H3 


2.8 ±0.8 


1.7 


2.5 


1.7 







1 See Sect. |5.l| for the e xpla nation of model PI, Sect, 
for P2 & P3, and Sect. O for El & E2 



5.3 



Figure ^ shows the cumulative contributions of clus- 
ters, type 1 AGNs and type 2 AGNs in model PI to the 
H / sh(0)(/ s }(/h) values for selected correlations. The ob- 
served value with ±1<j range is shown for each correlation. 
As we see in Table || and Fig. [| model PI is consistent 
within 2cr with the observed correlations between A2 maps 
and bright R6 sources (T3 and H3). We also see that these 
correlations are dominated by clusters of galaxies. On the 
other hand, the observed correlations between the bright 
source removed R6 map and the A2 maps are about a 
factor of 3 - 4 larger than the model predictions. In any 
case, the correlation is dominated by sources in the local 
universe (z 55 0.2) and the expected correlation does not 
depend on the cosmological evolution significantly. The 
contribution of galactic stars, effects of minor modifica- 
tions to the model, and the effects of possible extended 
components are discussed in the followings. 

5.2. Contribution of Galactic Stars 

Among the 835 serendipitously found X-ray sources in 780 
deg 2 of high galactic latitude sky in the Einstein Extended 
Medium Sensitivity Survey (EMSS) , about 25% are inden- 
tified with galactic stars (Stocke et al. 1991). The identifi- 
cation was 96% complete and the flux limits of the survey 
range between 0.5 - 34 x 10 _13 erg s _1 cm~ 2 in the 0.3-3.5 
keV band. When converted to the ROSAT R6 band, this 
corresponds to a flux range just below the detection limit 
in our source detection discussed in Sect. 3.2. Based on 
the flux-number (log N - log S) relation for the EMSS 
sample of stars, we have estimated the contribution of 
Galactic stars to the present correlation signals. Assum- 
ing a kT = 2 keV Raymond & Smith plasma spectrum 
with cosmic abundances, which is typical of a variable 
hard component of stellar coronal emission, the estimated 
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Fig. 3. The cumulative contributions of clusters, type 1 and 2 AGNs, and the total of these for model PI to the observed CCFs 
are plotted. Each component is drawn in a different line style as labeled in the lower left panel. The observed values with la 
error bars are also shown. The horizontal positioning of the observed data does not have any significance but placed there only 
for the convenience of the display. 



contribution to the R6 background versus A2 TOT corre- 
lation (T2) is w 3% and R6 source versus A2 TOT cor- 
relation (T3) is k 10% of the observed. Thus, the contri- 
butions of galactic stars to the observed cross-correlations 
arc within the statistical errors and not significant. 



5.3. Unabsorbed Soft Component associated with AGNs 

In our model PI, we used the AGN population synthe- 
sis model by Comastri et al. (1995). We here consider two 
minor modifications to this model. If a significant fraction 
of AGNs are embedded in groups of galaxies with X-ray 
emitting gas with kT ~1 keV, this could contribute to the 
soft-hard correlation. In many type 2 AGNs, there is an 
unabsorbed X-ray emitting component, which is believed 
to be electron-scattered X-rays from the nucleus. Comas- 
tri et al. (1995) et al. deliberately neglected this compo- 
nent because its luminosity is much weaker than the soft 
X-ray luminosity in type 1 AGNs with a comparable op- 
tical luminosity. We see the impact of these effects to our 
correlation. 



In model P2, each AGN has an unabsorbed Raymond 
& Smith plasma component of kT = 1 keV with a 0.5 - 
2 keV luminosity of 2 x 10 42 /i^ 2 ergs _1 . This is a typical 
value for most luminous groups of galaxies found in the 
ROSAT North Ecliptic Pole Survey (Henry et al. 1995). 
The spatial number density of groups of galaxies at this lu- 
minosity is several times 10~ 5 h§ Mpc~ 3 , which is compa- 
rable to the number density of Seyfert galaxies and about 
two orders of magnitude smaller than the number density 
of bright spiral galaxies. Thus, it is unlikely that every 
AGN is associated with a group of galaxies at this lumi- 
nosity. The purpose of this model is not to be realistic but 
to examine the effect of the additional component. Table |^ 
shows that adding this component has little effect on the 
correlations T2 and H2, thus fails to explain the excess 
correlation even in this extreme case. 

As another model, P3, we put a soft component to 
AGNs with fluxes proportional to the intrinsic (i.e. un- 
absorbed) luminosity. The underlying idea of this model, 
in terms of the unified scheme of AGNs, is to add an 
unbasorbed X-ray component in Seyfert 2, which is gen- 
erally considered as electron-scattered nuclear emission. 
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Like Madau et al. (1994) considered in their model, we as- 
sume that 2% of the nuclear emission in type 2 AGNs are 
from the scattered component and can reach us without 
being absorbed. As we see in Table |], this component is 
not sufficient to explain the observed excess correlation 
over model PI. 

5.4- Low-Luminosity Flat X-ray Spectrum Galaxies 

The Comastri et al. (1995) model deliberately neglected 
the contribution of low X-ray luminosity objects (L x Si 
10 42 /i^ 2 ergs- 1 ) to the hard X-rays (E > 2 keV), while 
they left some room for these objects to contribute to 
soft X-rays (E < 2 keV). However, Yaqoob et al. (1995) 
recently reported that NGC 3628, a low X-ray luminos- 
ity object, which had usually been classified as a star- 
burst, has a rather flat X-ray spectrum (a photon in- 
dex of a p h = 1.2). Based on this observation, they ar- 
gued that a class of sources represented by this object 
could contribute significantly to the XRB. Here we esti- 
mate the contribution of the possible low luminosity pop- 
ulation to our correlations. For an estimation, we con- 
sider low-luminosity X-ray galaxies of the 2-10 keV lu- 
minosity range of 10 40 — 10 42 /i^~ 2 ergs _1 at the present 
epoch. We assumed a power-law differential lumninos- 
ity function (d<i>(L x )/dL x oc £~ 7 ) with 7 = 1 and to- 
tal 2-10 keV present-epoch volume emissivity of 2 x 
10 38 /i5oergs _1 Mpc~ 3 . This volume emissivity is an upper 
limit from the low-luminosity source contribution derived 
by Miyaji et al. (1994) based on the XRB- IRAS galaxy 
cross-correlation and has been considered in the Yaqoob et 
al. baseline model. For this population, we have assumed a 
power-law spectrum with a p h = 1.2 and simple luminosity 
evolution of oc (1 + z) 3 - 2 up to z = 2. The contribution of 
this population to our correlation calculated using above 
assumptions is only ~ 1% for either T2 or T3. Note that, 
in this model, most contribution to the correlation signal 
comes from around z ~ 0.1, unlike the contribution to the 
total surface brightness, which is mostly contributed from 
the high-redshift region given this evolution law. Thus this 
estimate is insensitive to the detailed behavior of the cos- 
mological evolution. The above estimation shows that the 
low X-ray luminosity galaxies cannot explain the observed 
excess correlation in T2 or H2. 

5.5. Extended Emission Component 

We consider here the effect of the clustering term (or ex- 
tended components) to the observed W s h(0). In the ex- 
tragalactic component of the R6 (1 keV) map, Soltan 
et al. (1995) found a significant extended component in 
the auto-correlation function (ACF) at a scale of 9 5s 5°. 
They argued that the ACF signal they have found should 
be dominated by real extended structures rather than clus- 
tering of sources, based on comparison with the ACF in 
Soltan and Hasinger (1994) at smaller scales. They spec- 



ulated that this component could be due to the clustering 
of groups of galaxies. However, information on statistical 
properties of groups of galaxies was very scarce and they 
did not make further arguments on this possibility. 

The expected R6 ACF at ~ 1° - 3° from the clusterings 
of AGN and clusters used in model PI using Eq. 12 (in this 
case, two bands in the expression are the same) are two or- 
ders of magnitude smaller than the actual values observed 
by Soltan et al. (1995). We also checked if the AGN and 
cluster clusterings can explain the excess correlation in 
our zero-lag CCFs (T2 and H2). The expected clustering 
effects to the ROSAT - HEAO-1 zero-lag CCFs (because 
of the beam smearing, the zero-lag values are mainly due 
to the clustering at a scale of 1°) for the same model are 
also two orders of magnitude smaller than the observed 
excesses. For these estimates, we have assumed a form 
£(r,z) = (r/r )- 7 (l + z)~ 3 - £ for r < 3r and €(r,z) = 
for r > 3ro. The parameters used for the estimations are 
(r [kms" 1 ], 7) = (1000, 1.8), (2400, 1.8) and (880, 2.2) for 
the AGN-AGN, cluster-cluster (Bahcall 1988) and AGN- 
cluster (Lilje & Efstathiou 1988) correlation functions re- 
spectively. At the angular scales of degrees, the expected 
correlation is insensitive to the cosmological evolution of 
clustering and we have used e = 0. Since the discrepan- 
cies are two orders of magnitude, we also expect that the 
clustering of low X-ray luminosity objects (galaxies, low- 
luminosity AGNs), which have not been considered in the 
above calculation, cannot make much contribution, since 
their 2-10 keV X-ray volume emissivity should be lower 
than that of AGNs (Miyaji et al. 1994). Thus we expect 
that the Soltan et al.'s extended ACF component truely 
represents a new population of X-ray structure. 

Here we consider a picture that the component re- 
sponsible for the Soltan et al.'s ACF component is also 
responsible for the excess zero-lag CCF in T2 and H2. 
We constrain the nature of this component by compar- 
ing the upper-limit to our 9 = 2?5 CCF with their ACF 
and also investigate the contribution of this component 
to our zero-lag CCF. Since the spatial resolution of their 
ACF measurement is much higher than our CCF measure- 
ment, we have smoothed their ACF with the beam of our 
CCF (see Eq. [11]) and compared it with our CCF. The 
smallest angle where there is no overlap of the beams is 
at 9 ~ 2.5°. At this angle, the smoothed ACF value is 
~ 1 x 10~ 3 , while the 2a upper-limit of our CCF (T2) is 
3 x 10 -4 . Assuming that the detected R6 ACF is domi- 
nated by a single component with a single spectrum, the 
hard-to-soft flux ratio Ih/L s of this component is: 

(J h /J B ) = (W sh (I h ))/(W ss {I s )). (13) 

The upper limit to this ratio is 1.3 (TOT cts s _1 ) (R6 cts 
s -1 ) -1 . This corresponds to a Raymond & Smith plasma 
with kT 2.5 keV or a power-law photon index of a p h <J 
2.3 for a galactic column density of Nn = 1.8 x 10 20 cm -2 . 

Assuming that the structure causing the extragalactic 
1 keV (R6) ACF in Soltan et al. (1995) in the (0.1° < 9 < 
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3°) comes from a single extended component, we have cal- 
culated the expected zero-lag R6-A2 CCF (W sh (0)). These 
values are shown under the models El (assuming a kT = 
2 keV Raymond-Smith spectrum) and E2 (kT = 1 keV) 
(E represents the extended). A heavy metal abundance of 
0.4 (cosmic) has been used in these estimations. Since the 
ROSATR6 (1 keV) band is dominated by the Fc L emis- 
sion lines, this estimate is sensitive to the heavy element 
abundance and should be taken as a rough measure. 
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Fig. 4. The observed angular CCFs for the source-removed 
cases are shown with two models as indicated. The la error 
bars are from the scatters in the 'rotated' correlations (dotted 
lines in FigJ^). Note that the errors in data points are not 
independent but highly correlated 

Figure |4| shows the comparisons of the observed and 
modeled angular CCF for the source-removed cases. The 
models shown are El and E2 in addition to PI as la- 
beled. As we see in Fig. |], model (Pl+El) could explain 



the zero-lag amplitude of the CCFs. This means that if 
the extended component in the R6 ACF is a thermal 
plasma with kT ~ 2 keV this would also explain the ex- 
cess zero-lag CCF over the model PI prediction in our 
source-removed CCFs within 2a. Considering that the H2 
zero-lag correlation still exceeds the model Pl+El, while 
statistical significance is poor, a flatter spectrum in the 
E 2 keV range for this component is preferred. It is 
possible that the R6 ACF signal Soltan et al. (1995) found 
is composed of more than one component. For example, 
a composite model with a small-scale (6 Ss 2°) hard com- 
ponent and a large-scale (2° Ss 9 < 5°) soft component 
would better explain the data. 

6. Concluding Discussions 

As we see from the analysis presented above, AGNs 
and clusters account well for the correlation between the 
ROSAT bright R6 sources (^ 0.08 R6 cts s _1 , corresponds 
to 2.5 x lO'^ergs^ 1 cm~ 2 ) and the HEAO-1 A2 surface 
brightness. This correlation is dominated by the Pois- 
son effect from individual clusters of galaxies in the lo- 
cal universe. On the other hand, the correlations between 
the residual R6 background and the HEAO-1 1 A2 sur- 
face brightnesses are about a factor of 3 - 4 larger than 
the Poisson effect expected from the population synthe- 
sis model by Comastri et al. and clusters of galaxies based 
on the known X-ray luminosity functions. We have con- 
sidered a few possible components which have not been 
considered in the population synthesis model, i.e., AGN- 
group association, unabsorbed (scattered) soft component 
in type 2 AGNs, and low-luminosity X-ray galaxies. Any 
of these three are not likely to explain the observed excess 
correlation. 

Due to the large collimator beam of the HEAO-1 A2 
experiment, the zero-lag correlation could also be con- 
tributed by the clustering effect or an extended com- 
ponent. The extragalactic 1 keV (R6) angular auto- 
correlation function (ACF) by Soltan et al. (1995) clearly 
shows an extended structure up to a separation of sev- 
eral degrees. The origin of this component is not clear at 
this moment. They showed that the effects of rich clusters 
including the cluster-cluster clustering and the putative 
~ 10 Mpc-scale diffuse gas component in their model to 
explain the Abell cluster - 1 keV cross-correlation function 
only amounts to about one third of the R6 ACF signal. 
They suggested diffuse emission from groups of galaxies 
could possibly contribute to the R6 ACF. However, sta- 
tistical properties concerning poor groups, e.g. luminos- 
ity/temperature functions and clustering properties with 
rich clusters/AGNs, are scarcely known. Comparing the 1 
keV ACF with the upper-limit to our R6 - A2 (TOT) CCF 
at 9 = 2.5°, the minimum angular separation where beams 
do not overlap, we have constrained the spectrum of the 
extended structure causing this R6 ACF signal. This cor- 
responds to a Raymond & Smith plasma of kT <, 2.5 keV 
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or a power-law photon index of a p h il 2.3. Suppose the 
component responsible for the 1 keV ACF has a spectrum 
with an 'equivalent color-temperature' of kT ~ 2 keV, this 
component explains the correlation in excess of the pop- 
ulation synthesis model expectation, assuming that the 
extended R6 ACF structure is composed of a single com- 
ponent. A picture where there are a harder small-scale 
component and a softer large-scale component would bet- 
ter explain the data. 

There remains also a possibility that some major 
change to the Comastri et al. population synthesis model 
or an introduction of a new population is needed, while 
the extended component has a lower temperature (or a 
steeper spectrum). The current observation cannot dis- 
criminate between these two, since it is limited by the 
low spatial resolution of the HEAO-1 A2 experiment. A 
large area survey at E ^ 2 keV with higher spatial reso- 
lution in a future mission such as ABRIXAS (Friedrich et 
al. 1996) would enable us to discriminate between these 
two pictures. With an analysis similar to that presented 
here, but with a higher spatial resolution, we can con- 
strain the spectrum of the extended component and give 
another crucial observational constraint to the population 
synthesis models simultaneously. 

Acknowledgements. TM is supported by a fellowship from the 
Max-Planck Society. TM appreciates the hospitality of the As- 
trophysikalisches Institut Potsdam during his visit. The au- 
thors thank Keith Jahoda for his help in handling the HEAO- 
1 A2 data and Elihu Boldt for encouragements. This work 
has made use of data from two X-ray satellites, ROSAT and 
HEAO-1 , in orbit about a decade apart from each other. The 
authors appreciate the effort of all people which lead to the suc- 
cess of these two missions. Thanks are also due to the referee, 
Luigi Danese, for useful comments and suggestions. 

References 

Allen J., Jahoda K.,& Whitlock, L. 1994, Legacy 5, 27 
Awaki H. Koyama K., Inoue H., & Halpern J.P. 1991, PASJ 
43, 195 

Bahcall N.A. 1988, ARA&A 26, 631 

Barcons X., Branduardi-Raymont C, Warwick, R.S., Fabian, 

A.C., Mason K. O., McHardy I., & Rowan- Robinson M., 

1994, MNRAS, 268, 833 
Barcons X., Franceschini A. Danese L., De Zotti G. & Miyaji 

T. 1995 ApJ, 455, 480 
Boyle B.J., Griffiths R.E., Shanks T., Stewart G.C. & Geor- 

gantopoulos I. 1993, MNRAS 260, 49 
Boyle B.J., Shanks T., Georgantopoulos I., Stewart, G.C, & 

Griffiths, R.E. 1994, MNRAS 271, 639 
Carrera F.J. et al. 1991, MNRAS 249, 698 
Comastri A., Setti G., Zamorani G., & Hasinger G. 1995, A&A 

296, 1 

Danese L., Toffolatti L., Franceschini A., Martin-Mirones J.M., 

& De Zotti G. 1993, ApJ, 412, 56 
David L. P., Slyz A., Jones C, Forman W., Vrtilek S. D., & 

Arnaud, K. A. 1993, ApJ 412, 479 



11 

Edge A. C, Stewart G. C, Fabian A. C, & Arnaud K. A. 

1990, MNRAS 245, 559 
Friedrich P. et al. 1996, in Rontgenstrahlung from the Universe, 

eds. Zimmermann H.U. & Triimper J., MPE report, in press 
Gioia I.M., Maccacaro T., Schild R.E., Wolter A., Stocke J.T., 

Morris S.L. & Henry J.P. 1990, ApJS, 72, 567 
Hasinger G., Burg R. Giacconi R. Hartner G., Schmidt M., 

Triimper J., & Zamorani G. 1993, A&A 275,1 
Henry J.P. et al. 1992, ApJ 386, 408 
Henry J.P. et al. 1995, ApJ 449, 422 
Jahoda K. 1993, Adv. Space Res. 13, No. 12, 231 
Jahoda K., Lahav O., Mushotzky R., Boldt E. 1991, ApJ 378, 

L37 

— . 1992, ApJ 399, L107 (Erratum) 
Lahav O. et al. 1993, Nature 364, 693 
Lilje P.J. & Efstathiou G. 1988, MNRAS 231, 635 
Madau P., Ghisellini G. & Fabian A.C. 1994, MNRAS 270, L17 
Miyaji T. 1994, PhD Thesis, University of Maryland 
Miyaji T., Lahav O., Jahoda K., & Boldt E., 1994, ApJ 434, 
424 

Mulchaey J.S., Mushotzky, R.F., & Weaver, K.A. 1992 ApJ 
390, L69 

Pfeffermann, E. et al. 1986, Proc SPIE 733, 519 
Piccinotti G. et al. 1982, ApJ 253, 485 

Roche N., Shanks, T., Geaorgantopuolos, I., Boyle B.J., & Grif- 
fiths, R.E. 1995, MNRAS 273, L15 

Rothschild R.E. et al. 1979, Sp. Sci. Instr., 4, 269 

Schartel N, Schmidt M., Fink H. H., Hasinger G., & Triimper 
J. 1995, preprint 

Shafer R. A. 1983 PhD Thesis, University of Maryland 

Soltan A. & Hasinger G. 1994, A&A 288, 77 

Soltan A., Hasinger G., Egger R., Snowden S., & Triimper J. 
1995, A&A, in press 

Snowden S.L. & Schmitt, J.H.M.M. 1990, Ap&SS, 171, 207 

Snowden S.L., McCammon D., Burrows D.N., Mendenhall J. A. 
1994, ApJ 424, 714 

Snowden S.L. et al. 1995, ApJ in press 

Stocke J.T., Morris S.L., Gioia I.M., Maccacaro T., Schild R., 
Wolter A., Fleming T.A., & Henry, J.P. 1991, ApJS 76, 813 
Triimper J. Adv. Space Res. 2, No.4, 241 

Voges W. 1992, in Proceedings of Satellite Symposium 3, Space 
Science with particular emphasis on High-Energy Astro- 
physics, ed. T.D. Guyenne & J.J. Hunt (noordwijk:ESA 
Publication Division), 9 

Voges W. 1995, private communication 

Yaqoob T., Serlemitsos P.J., Ptak, A., Mushotzky R., Kunieda 

H., & Terashima Y. 1995, ApJ, 455, 508 
Zimmermann H.U., Becker W., Belloni T., Dobereiner S., Izzo, 

C, Kahabka P. & Schwentker O. 1994, EXSAS User's 

Guide, MPE Report 257, Garching 



This article was processed by the author using Springer- Verlag 
DT E X A&A style file L-AA version 3. 



